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OBLIQUITY HISTORIES OF EARTH AND MARS: INFLUENCE OF INERTIAL AND 
DISSIPATIVE CORE-MANTLE COUPLING; Bruce G. Bills, LPI, Houston, TX 

For both the Earth and Mars, secular variations in the angular separation of the spin axis 
from the orbit normal are suspected of driving major climatic changes (1, 2,3,4). There is thus 
considerable interest in determining, as accurately as possible, the amplitude and timing of 
these obliquity variations. If the orientation of the orbital plane were inertially fixed, and the 
planet were to act as a rigid body in its response to precessional torques, the spin axis would 
simply precess around the orbit normal at a fixed obliquity e and at a uniform angular rate a 
cos(e). The precession rate parameter 

a = 3 (C-(A+B)/2) E /V(bj) 3 PI 

2 Cw 

depends on the principal moments of inertia (A<B<C) and rotation rate w of Ihe pei rt urbed 
body, and on the gravitational masses /i = GM and semiminor axes b = a(l-e ) o te 
perturbing bodies. For Mars, the precession rate is not well known, but probably lies in the 
interval 8-10 arcsec/year (5,6,7,8,9,10). The much larger precession rate for the Earth (-50.3 
arcsec/year) is partly due to greater proximity to the Sun (which contributes roughly 1/3 of the 

total), but is mostly due to the presence of the Moon (11). 

Gravitational interactions between the planets lead to secular motions of the orbit planes 

of the form 


sin(I) e in = E Nje^j 1 + 5 j> 


[ 2 ] 


where I and fi are the inclination and longitude of the node, and Nj, sj and $j are amplitudes, 
rates and phase constants (12,13). In the rigid body case, the spin axis still attempts to precess 
about the instantaneous orbit normal, but now the obliquity varies. A first order solution for 
the obliquity can be written in a form similar to [2] but with amplitudes KjNj, where the 
admittance Kj has the value (14,15,16) 


Kj = 




sj + a cos(e) 


[3] 


As the orbital precession rate constants s: are all negative and fall within the range {-26.3 
arcsec/year < s: < 0), while the spin precession rate constants are positive, the potential exists 
for significant resonant amplification of the obliquity if the denominator of [3] approaches zero 
(17,18). The possibility of actual singularities in [3] is not too worrisome, as the linear analysis 
leading to that form of the admittance is no longer applicable in the immediate vicinity of a 
resonance. However, the physical model of rigid rotation which leads to that formula is almost 

certainly too simplistic. _ . . . . 

The hydrostatic figure of a planet represents a compromise between gravitation, which 
attempts to attain spherical symmetry, and rotation, which prefers cylindrical symmetry (19). 
Due to their higher mean densities, the cores of the Earth and Mars will be more nearly 
spherical than the outer layers of these planets. The direct gravitational torques on the core will 
thus be inadequate to make it precess at the same rate as the mantle. For the Earth, where the 
structure is relatively well known, the core oblateness is only about 3/4 that required for 
coprecession with the mantle (20). However, it is clearly the case that the core and mantle 
precess at very nearly the same rate (21,22). Two different types of torques contribute to the 
coupling. 
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On short time scales it is appropriate to consider the core to be an invisid fluid 
constrained to move within the ellipsoidal region bounded by the rigid mantle (23,24). The 
inertial coupling provided by this mechanism is effective whenever the ellipticity of the 
container exceeds the ratio of precessional to rotational rates. If the mantle were actually rigid, 
or even elastic, this would be an extremely effective type of coupling. However, on sufficiently 
long time scales, the mantle will deform viscously and can accomodate the motions of the core 
fluid. The inertial coupling torque exerted by the core on the mantle will have the form Tj = 
^i [ w m x w cJ- A fundamentally different type of coupling is provided by electromagnetic or 
viscous torques (25,26). The dissipative coupling torque exerted by the core on the mantle will 
have the form T d = -k d [w m - w c ], This type of coupling is likely to be most important on 
longer time scales. In each case, the mantle exerts an equal and opposite torque on the core. 

The admittance which relates inclination amplitude to obliquity amplitude is now a 
complex quantity which can be written in the following form (27,28) 

K(s) = sfs ± a m ± r?) [4] 

(s + a m )(s + a m + ri) - Aa i 7 m 

The coupling constants are rescaled (£ m = k d /C m , 7 m = w cos(e) kj/C m ) and then combined 
to form a single complex parameter t? m = 7 m - i/V Also, a m = “m cos(e) + cos(I) dfi/dt. Core 
parameters (a c , 0 c , 7 c ,i? c ) are defined analogously, and the sum of core and mantle parameters ri 
= + *lc * s unsubscripted, and Aa = a m -a c . Viscous relaxation of the mantle is included 

by multiplying both 7 values by (isr)/(l+isr), where r is the effective Maxwell relaxation time 
of the mantle. 

There are several features to note, in comparing [3] and [4], As is typical of a forced 
oscillation with damping, the response now lags the forcing by an amount which depends on 
the frequency of forcing and the strength of the viscous coupling. The most evident difference 
in response is near resonance, since [4] exhibits no singularities. However, even away from 
resonance, the inclusion of possible differential precession can modify the obliquity history by 
amounts that could have climatic significance. Unfortunately, the coupling constants are not 
well known, even for the Earth, and are almost completely unconstrained for Mars. Thus, some 
caution is advised in constructing climatic history scenarios which depend on details of the 
obliquity history. 
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Implications of the low carbonate abundance in the optical surface of Mars. Diana L. Blaney. 

Planetary Geosciences Division, Hawaii Institute of Geophysics, 2525 Correa Rd. Honolulu, HI 
96822. 

There are two major reasons for expecting carbonates on Mars: 1) if Mars underwent a period of climatic wanning 
due to a thicker C02 atmosphere, carbonates would have formed at a rapid rate (1,2) and 2) the thermodynamically 
stable weathering products of igneous material on Mars under current conditions also include carbonates (3). 
Earth-based spectroscopy of the Martian surface in the 4.0 pm region has placed an upper limit on the carbonate 
abundance of < 1 - 3 wt% in the optical surface of the regions measured based on laboratory mixing models (4,5). 
The implications for weathering and climate change on Mars using various models to explain the low carbonate 
abundance are presented. 

Model 1. The observed upper limit on carbonate abundance (he. <1 wt%) is representative 

of the Martian regolith. A 1 km mega-regolith with 1 wt% carbonate could at most contain 0.6 bars of C02 

(2), well short of the 3-5 bars (2,6,7) needed for global wanning by a C02 greenhouse at average orbital parameters. 

A consequence of an early wet warm Mars is that liquid water would be present to remove C02 from the atmosphere 
by formation of carbonates would occur at rates of approximately 1 bar C02 removed every 10 7 years (1,2). If the 
low carbonate abundance observed (0.6 bars), accurately reflects the planetwide inventory, then a 3-5 bar early 
greenhouse lost most of its C02 by other mechanisms and had a duration of a few million years. A short greenhouse 
contradicts the geomorphic evidence which shows valley networks (a major geomorphic feature supporting an early 
greenhouse), forming during throughout the first billion years. Recycling carbonate to replenish the atmospheric 
C02, whether by volcanic (2) or by impact (8) processes, may be able to extend the length of the greenhouse, but the 
last few bars would have to have been stored in the near surface environment (i.e. if you are forming carbonates by 
surface water, the last carbonates formed must be near the surface). 

A low carbonate abundance would also imply that the weathering products observed (i.e. the dust) are not the 
thermodynamically stable end-members of weathering under current ambient Martian conditions. As an illustration, 
consider the thermodynamically stable decomposition products from the weathering of basalt of the same composition 
as Shergotty. Shergotty is a basaltic achondrite from the SNC (S=Shergotty, N=Nakhla, and C=Chassigny) group of 
meteorites. SNC are thought for a variety of reasons to come from Mars (e.g. 9). The weathering products of a 
Shergotty-like basalt which has reached thermodynamic equilibrium by gas-solid interactions at 240°K and current Mars 
atmospheric conditions, based on (3) indicate a total carbonate abundance of 32.48 wt %. If the 1 wt% carbonate in 
the dust is a reasonable upper limit, then about only about 3 % of the dust could have formed by gas solid reactions 
carried to equilibrium. The majority of the dust would have formed by other processes. Similarly, if liquid water 
containing dissolved C02 and 02 from the current Mars atmosphere were involved, again based on (3), the end 
products would contain 18.07 wt% carbonate, leaving 5.5 wt% of the dust to form by this process. Clearly, either the 
majority of the dust formed by processes other than surface weathering under current ambient conditions, or the 
formation process has been halted mid-way by kinetic barriers leaving metastable products. 

Model 2. Carbonate deposits are located at depth. Carbonates deposits at depth are impossible to 
detect with spectroscopy, as spectroscopy looks only at the optical surface. Previous models calling for the recycling 
of carbonates to sustain an early greenhouse (2,8) assume that carbonates are well mixed in the regolith. Carbonate 
recycling rates from both impact and volcanic burial are too slow to replenish the carbon dioxide fixed as carbonates 
by liquid water (8). Carr (1989) suggested that the rate of carbonate fixation may have been overestimated. If liquid 
water was not continuously present on the surface, then carbonate formation would have preceded at a much slower rate. 
Temperatures, while warmer than the current epoch, would have remained just below the point at which liquid water 
could exist. Ground water would be close to the surface and under proper hydrostatic conditions would break out and 
form valley networks. Even with the slower carbon dioxide fixation rates, carbonates would form in the near surface 
environment where they should be detectable. Thus, a mechanism to for locating carbonates at depth is needed. 

An alternative hypothesis calls for stratified carbonate layers. If carbonates exist at depth in stratigraphic ally 
localized deposits a possible formation scenario could be: a) Mars is warm, valley networks form; b) carbonates form 
rapidly, decreasing the atmospheric C02 pressure until temperatures drop below the temperature needed to keep liquid 
water around to form carbonates; c) carbonate formation ceases and volcanic and impact processes bury the carbonates; 
d) carbonates are buried to a depth where rapid decrepitation occurs, releasing a pulse of C02; e) the released C02 
causes another episode of greenhouse warming and valley network formation, and carbonates begin to form again. If a 
carbonate decrepitation depth of about 17 km, as suggested by Carr (1989), is valid then the cycle of rapid carbonate 
formation followed by burial could repeat itself about 10 times in the first billion years according to liquid stirring 
models (8). Declining impact rates would then slow the burial rates, leaving the carbonate layer stratified somewhere 
in the Martian crust. There is some geomorphic evidence for a strength discontinuity at about 1 km depth (10) which 
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could be a carbonate layer. Geochemical arguments which focus on the calcium deficiency of the Viking fines when 
compared to SNC meteorites also indicate that a carbonate layer may be present at depth (11). 

Model 3. The low carbonate abundance is the result of secondary weathering processes 
which have removed carbonates from the regoiith. The presence of the sulfate in the Martian fines 
measured by Viking Landers 1 and 2 (12), may have major implications for carbonate abundance. 

The stability of carbonates in the presence of gaseous S02 has been addressed by (13) and (14) with differing 
results. Clark et. al. (13) performed laboratory experiments which indicated that carbonates decompose readily in the 
presence of S02 leaving sulfates. Sidorov and Zoltov (14) examined kinetic constraints and showed that the reaction 
encounters a kinetic barrier which would make this reaction negligible. Understanding the possible role of S02 in 
destroying carbonates is critical in understanding the composition and evolution of the fines at the Martian surface. 
Alteration of the surface by airborne reactive components may be a major process in the formation of salts on Mars. 

Settle (15) proposed a sulfate aerosol model for the formation of the sulfate duracrust at the Viking sites. The 

acidic nature of these aerosols would tend to decompose carbonates readily. If carbonates do decompose in the presence 
of S02 vapor or sulfate aerosols, C02 would be released from the decomposed carbonates back into the Martian 
atmosphere. A terrestrial example of aerosol formed salts are the sulfate and nitrate deposits form in Antarctica from 
the deposition of stratospheric anions, which react with local rocks which provide the cations to form these deposits 
(16). 

Carbonates and sulfates must be considered together. If carbonates are absent from the surface due to interactions 

with volcanic aerosols, then the rate of S02 released by volcanoes on Mars must be greater than the rate of regoiith 

overturn exposing new carbonates. Fresh craters which have occurred after the height of volcanism would be expected 
to show evidence of carbonates brought up from depth. The rates of both volcanism and cratering were higher during 
the early portion of Martian history than at the present (e.g. 17, 18). If the sulfur in the Martian regoiith is the result 
of sulfate aerosol deposition throughout Martian history, then the release of C02 by sulfate aerosols may have played 
a critical role in the rapid recycling of C02 throughout Martian history. 

Model 4. Weathering environments on Mars did not form carbonates. The majority of materials 
examined on the Martian surface are weathering products formed by unknown mechanisms. As discussed previously, 
the expected carbonates from surface atmospheric interactions do not seem to be present in large quantities in the 
optical surface of Mars. Processes unique to Mars may have formed the observed dust. A weathering mechanisms 
which may have created the dust is impact ejecta from volatile rich targets (19). 

In addition to the hydrothermal alteration of material in throughout sheets, Kieffer and Simonds (19) proposed that 
a fine ash deposit formed by the rapid expansion of silicate melt/vapor by superheated steam from pore water or by 
carbon dioxide from decrepitated carbonates. The ash formed by the interaction of melt and vapor and would occur at 
the top of the crater’s stratigraphic sequence. Erosion has removed this ash layer from terrestrial craters which, by 
analogy with melt sheets from impacts into crystalline material, would be 1-5% of the volume excavated. Impacts 
into volatile rich material provide a way of producing a large volume of fine grained material of unknown composition. 
The alteration of the melt in the vapor cloud would be controlled by the dominant volatile propelling the expansion of 
the melt-vapor cloud. Thus, if water was the principle volatile, oxidation and hydration products would dominate and 
little if any carbonate would form. Due to the speed with which this event would occur, the altered material would be 
poorly crystalline or crystalline on a small scale. 
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Earth-Based Telescopic Observations of Mars in the 4.4 pm to 5.1 pm Region. Diana L. 

Blaney and Thomas B. McCord. Planetary Geosciences Division, Hawaii Institute of 

Geophysics, 2525 Correa Rd. Honolulu, HI 96822. 

Measurements. The telescopic measurements of Mars reported here were made on August 19, 1988 UT at the NASA 
IRTF facility. Spectra were obtained with the Cooled Grating Array Spectrometer (CGAS) which utilizes a thirty-two 
element InSb line array with a resolution of R=300. Measurements were made by taking data at two different grating 
positions, with an eleven channel overlap between positions. Segment 1 covered 4.405 pm - 4.86 pm, and segment 
2 covered 4.67 pm - 5.13 pm. Spectra were normalized to unity at 4.71 pm and were reduced using the star BS437 as a 
standard. A 2.7 arc-second aperture was used, producing a spot size 900 km in diameter at the sub-earth point. The 
measurements were taken as part of a sequence of measurements which also covered the wavelength region between 3.2 

pm and 4.2 pm using three grating positions discussed(l). 

The Martian atmosphere has a large number of atmospheric absorptions from both CO and C02 in the wavelength 
region between 4.4 pm and 5.1 pm in addition to significant thermal components. These complications must be 
considered in the interpretation of the spectra. 

Detection of an Absorption Feature at ~4.5 pm --Indication of Sulfates. A modeled atmospheric 
spectra (PC02 =7 mbar, T=240K, C0/C02=0.002, Airmass = 2.0) based on atmospheric models by Encrenaz is plotted 
over the spectra from the Argyre region (Figure 1). Note that all the atmospheric absorption features (shown with 
arrows) are clearly present in the telescopic data and that most features appear to be weaker than predicted by the 
model. The exception to this is the rise out of the 4.2 pm - 4.4 pm C02 absorption which has a much shallower slope 
in the telescopic data than in the model atmosphere spectra. In addition to this gradual slope, an inflection at 4.5 pm 
is present in all the spectra and is extremely well defined in the Eastern Solis Planum and Argyre spectra (Figure 2). 
These measurements indicate that a surface absorption, centered near 4.5 pm, is present on the Martian surface. 

The wavelength region between 4.4 pm and 5.1 pm contains sulfate absorptions, shown in figure 3a, for 
anhydrite, gypsum, and MgS04, with the the 2v3 vibrational overtone of the S04-2 anion being centered at 4.5 pm. 

The multiband structure seen in the sulfate spectra is produced by the metal cation complexing with the S04 -2 
lowering the symmetry of the sulfate. The lower symmetry causes degeneracy in the v3 and v4 modes, and the 
appearance of the vl and v2 modes which are not infrared active under the tetrahedral symmetry that the S04 -2 anion 
occupies. Unfortunately the Mars atmosphere cuts off the short wavelength portion of the absorption. The sulfates 
shown in figure 3a are shown in figure 3b at the resolution and wavelength of the telescopic data. Note that there does 
not appear to be any of the structure associated with the sulfates shown above in figures 3 a and b present in the 
telescopic data. 

The location of this feature at 4.5 pm, the 2v3 S04 overtone, is especially exciting as the Viking lander sites 
detected large quantities of sulfur (~7 wt% S03) in the Martian soils. The sulfur content was highly variable even at a 
local level. Sulfur abundance ranged from 5.9 wt% to 9 wt% S03 at Chryse and from 7.6 wt% to 8.5 wt% S03 at 
Utopia (2). 

Spatial Variation in the 4.5 pm band depth — Implications for sulfate distribution. In order to 
isolate atmospheric effects from surface absorptions we compared spectra which had the same 3.81 pm band depth. The 
3.81 pm band is caused by C02 with 160 and 170 (1). The isotopicaUy heavy C02 absorption provides an indicator, 
independent of possible sulfate absorption, of the similarity of atmospheric conditions of the various spectra. The 
spectra for the regions Valles Marineris, Argyre Basin, Eastern Solis Planum and Eastern Tharsis (figure 2) all have 
3.81 pm band depth between 3.1% and 3.3%. 

In comparing the Valles Marineris, Argyre Basin, Eastern Solis Planum and Eastern Tharsis spectra, the 4.5 pm 
inflection discussed above is probably the strongest indicator of sulfate abundance. In order of strength of deepest 
absorption to weakest the regions are ordered 1. Eastern Solis Planum, 2. Argyre Basin, 3. Eastern Tharsis, and 4. 
Valles Marineris, although in no spectrum does the rise out of the carbon dioxide band match the atmospheric model. 
A note of caution must be offered in taking this approach to estimating relative sulfate abundance as the strength of 
the 4.5 pm band could indicate not only the abundance but changes in mineralogy and degree of crystallization. 
However, the four regions do follow a progression due to either sulfate abundance or chemistry and there is significant 
variation between regions shown in figure 2. 

Conclusions. The rise out of the 4.2 pm - 4.4 pm C02 band cannot be matched by solely atmospheric constituents. 
A surface absorption must be added at roughly 4.5 pm in order to decrease the reflectance rise and produce die 4.5 pm 
infle ction which is present in the data. The known presence of sulfates on the Martian surface and the location of this 
feature at the 2v3 overtone of the S04 -2 anion indicate that the surface absoiption is probably caused by sulfates on 
the Martian surface. An exact match to a terrestrial sulfate mineral has not been made but it is suggested that the 
mineral has very weak structure and thus a high degree of symmetry . Significant variation exists between the spectra 
at 4.5 pm. In order of strength of deepest absorption to weakest the regions are ordered Eastern Solis Planum, Argyre 
Basin, Eastern Tharsis, and Valles Marineris for the four regions measured at similar Mars atmospheric conditions. 
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SPECTROSCOPY ON Mars — SlfLEAXES 





^Figure 2. Spectra for four regions with similar 
-J 3 . 8 I pm band depths. Regions are Eastern Solis 
] Planum, Ridged Plains, Heavily Cratered Terrain; 
J Axgyre Basin; Eastern Tharsis; and Valles 
j Marineris. Spectra are arranged in order of 4.5 
l pm band depth. 
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Figure 3 a. Reagent MgS04, gypsum, and 
anhydrite reflectance spectra. Spectra are relative 
to sulfur and have a grain size of <34 pm. 
Spectra measured at 2 cm "1 resolution, b. 
Reagent MgS04, gypsum, and anhydrite 
reflectance spectra. Spectra are relative to sulfur 
and have a grain size of <34 pm. Spectra 
convolved and interpolated to wavelengths of the 
telescopic observations. 
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RADIATION DAMAGE IN GERMANIUM DETECTORS: IMPLICATIONS FOR THE 
GAMMA-RAY SPECTROMETER OF THE MARS OBSERVER SPACECRAFT; J. Bruckner, 
M. Korfer, and H. Wanke, Max-Planck-Institut fur Chemie, Mainz, F. R. Germany; A. N. F. 
Schroeder, Universitat zu Koln; F. R. Germany; D. Filges and P. Dragovitsch, Institut fur 
Kemphysik, KFA Jiilich, F. R. Germany; P. A. J. Englert, San Jose State University, CA, USA; 
R. Starr and J. I. Trombka, NASA Goddard Space Flight Center, MD, USA; I. Taylor, Princeton 
Gamma-Tech, NJ, USA; D. Drake and E. Shunk, Los Alamos National Laboratory, NM, USA. 

High-purity germanium (HpGe) detectors will be used in space missions to measure gamma- 
rays emitted by planetary bodies or astronomical sources. The upcoming NASA Mars Observer 
mission, scheduled for launch in 1992, will explore the surface of the planet Mars for 1 Martian 
year. One of the scientific instruments on-board the spacecraft will be a gamma-ray spectrometer 
(GRS) that will utilize a coaxial n-type HpGe detector. The GRS shall make high-resolution 
energy measurements of gamma rays emitted by the planetary surface. Cosmic-ray bombardment 
of the planetary surface induces gamma-ray emission due to subsequent primary and secondary 
nuclear interactions. Gamma-ray measurements made aboard the polar orbiting Mars Observer 
spacecraft can be used to infer chemical composition of the Martian surface. 

Previous experiments [1] and flight experience [2] have indicated that HpGe detectors 
exposed to cosmic radiation can suffer significant radiation damage. As a result, the energy resolu- 
tion of the HpGe detector will degrade to such a degree that it is not useful for further measure- 
ments. 

The purpose of this experiment was to obtain more detailed information on the behavior of 
HpGe detectors under irradiation with high-energy charged particles in space and to test methods 
of annealing radiation-damaged HpGe detectors in a closed system. Basic scientific as well as 
engineering data on detector performance under simulated space exposure and operating conditions 
were collected. These data will be used to help in the design, control, and understanding of how to 
maintain the health of HpGe detectors during future space missions. 

Several HpGe detectors were exposed to a proton irradiation at the French accelerator facility 
Satume (Laboratoire National Satume, Saclay). One detector was a flight-type version of the Mars 
Observer gamma-ray HpGe detector with a volume of 130 cm3. Six other HpGe detectors having 
volumes of 90 cm3 were grouped in pairs, each pair mounted in a special double-cryostat. Each 
pair consisted of a n-type and a p-type germanium crystal. 

Simulation calculations on the interaction of protons with germanium were earned out by 
using the HERMES code system. Implications from the calculations and space environment condi- 
tions led to special requirements for the proton beam and detector setup. At the accelerator a 1.5- 
GeV proton beam with a cross sectional area of 20 cm in diameter and with an intensity distribu- 
tion of less than 20 % variation was obtained. The proton flux was as low as about 105 protons 
cm- 2 sec-! and was measured with plastic scintillation counters and monitor foils. The HpGe detec- 
tors were arranged inside the beam cross section in such a way that no Ge crystal was in the geo- 
metrical shade of another. By using a special heating device the detectors were held at operating 
temperatures of 90 K, 100 K, and 120 K to cover a temperature range expected for a planetary 
mission. The irradiation of the detectors occurred in several steps until an accumulated fluence of 
about 108 protons enr 2 was reached. The detectors were irradiated in predetermined increments of 
proton fluence. During irradiation interruptions, gamma-ray spectra of 60Co were taken to evaluate 
the performance of the detectors as a function of accumulated charged particle fluence. 

At completion of the irradiation, all detectors were characterized and cooled stepwise to LN 2 
temperatures. Then, the detectors were transported in special containers to the laboratory at Mainz, 
where a detailed characterization of their performance was carried out 

A total proton fluence of 108 protons cm* 2 is about equivalent to one year of exposure in 
space. At elevated temperatures the observed energy resolution degradation after accumulation of 
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this fluence would not allow for high-resolution gamma-ray spectroscopy. But, as shown later, 
annealing can remove the damage. 

During the stepwise irradiation, the peak shape of the n-type detectors showed a significant 
change. Before the irradiation the shapes were either pure Gaussian or showed a little exponential 
tailing at the low-energy side. During the irradiation die exponential tailing evolved until it domi- 
nated the entire low-energy side of the the peak. A Gaussian shape was maintained at the high- 
energy side. More irradiation produced a new feature depending on the operating temperature: a 
broad Gaussian peak sitting on top of the exponential tailing (in a logarithmic display visible as a 
bump). 

The full-width-at-half-maximum (FWHM) of the three 90-cm3 n-type detectors is shown in 
Figure 1 as a function of proton fluence and operating temperature. The detector held at 90 K 
showed the best resistance against degradation of energy resolution. The detector held at 100 K 
(further on called "Yellow") exhibited an increase of degradation to 6 keV after being exposed to a 
fluence of 1.1x108 protons cm* 2 ; prior to the irradiation the resolution was 2.1 keV at 100 K. At 
an accumulated fluence of 6xl0 7 protons cm-2, which corresponds to about a half year of exposure 
in space, the energy resolution of detector Yellow was still below 3 keV. After the irradiation 
detector Yellow was cooled down to 90 K and its FWHM improved to 4.5 keV. The detector at 
120 K showed an almost linear dependence of FWHM vs. fluence in contrast to the curve of 
Yellow, which led to an resolution degradation from 2 keV before the irradiation to 3 keV at an 
accumulated fluence of 3xl0 7 protons cm-2. 

Temperature cycling of damaged 
detectors between a reference tempera- 
ture of 90 K and temperatures up to 
130 K proved that irreversible addi- 
tional loss in detector performance 
occurred above 120 K. For detector 
Yellow, it was observed that between 
125 and 130 K an irreversible increase 
of damage occurred, e. if the detector 
was cooled down to 90 K after being at 
130 K, its energy resolution had 
increased to 6.7 keV compared to 
4.5 keV at 90 K. Heating up to room 
temperature produced only a small fur- 
ther increase of the resolution degrada- 
tion (7.2 keV at 90 K). 

During annealing (heating up to 100°C or higher) the detectors stayed in their special 
designed cryostat (the Ge crystal was encapsulated in a vacuum-tight inner can). After having 
annealed detector Yellow at a temperature of 105°C for 11 hours, almost all damage was removed 
(resolution improved from 7.2 to 2.4 keV). Further annealing of 50 hours at the same temperature 
brought the energy resolution back to its original value. A further indication of complete recovery 
of the annealed detector was the lack of variation in the intrinsic detector energy resolution as a 
function of temperature in the range 90 K to 120 K. 

The irradiation and annealing cycles of the detectors simulated space exposure and operating 
procedures as closely as possible. The annealing temperatures applied to the detectors are available 
on the Mars Observer GRS. The study showed that HpGe detectors can successfully be used in 
long-term space experiments. 

References: [1] Pehl R.H., Madden N.W., Elliou J.H., Raudorf T.W., Trammell R.C., and Darken L.S. 
(1979) IEEE Trans. Nucl. Sci. NS-26, 321-323. [2] Mahoney W.A., Ling J.C., and Jacobson A.S. (1981) Nucl. 
Instr. Meth. 185, 449458. 
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CHEMICAL EVOLUTION AND OXIDATIVE WEATHERING OF MAGMATIC 
IRON SULFIDES ON MARS 

Roger G. Burns and Duncan S. Fisher, Department of Earth, Atmospheric and Planetary 
Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139 

INTRODUCTION. Ferric-bearing phases in the martian regolith testily to oxidative weathering reactions 
near the surface of Mars. A paragenetic link between primary magmatic sulfide minerals on Mars and their 
weathering products (gossans) has been proposed by us [1-3]. This model is based on petrogenetic 
associations between komatiitic rock-types. Viking geochemical data, SNC meteorites, and terrestrial 
magmatic sulfide ore deposits [4]. It appears that unique conditions during the evolution of the martian 
crust have conspired to generate an acidic groundwater (permafrost) environment that has facilitated 
chemical weathering of basaltic rocks on Mars. 

Some of the steps that have been proposed during oxidative weathering of iron sulfides [1-4] are 
summarized by reactions formulated in Table 1. Initially, deep-weathering reactions involving dissolved 
ferric iron in acidic groundwater percolating from the surface converts pyrrhotite to pyrite (reaction Hi). 
Pyrite is then oxidized either by ferric ions 121 or by groundwater equilibrated with oxygen in the 
atmosphere 121. even under an oxygen partial pressure of 10" 5 bar existing in the present-day martian 
atmosphere. Strongly acidic (pH 1 to 5) and sulfate-rich solutions are generated which stabilize and 
mobilize a variety of simple and complex ions of Fe, including Fe 2+ , Fe 3 *, FeS0 4 + , FeOH 2+ , etc. Above the 
water table, these dissolved ions become unstable and produce hydronium jarosite 141 or eventually 
goethite by further oxidation 1SL hydrolysis ifil, and decomposition 121 reactions, each of which 
liberates H + ions thereby promoting chemical weathering of basaltic silicate minerals 181- 

To test this weathering model, komatiitic pynhotites and olivines have been exposed to sulfuric acid 
solutions, with and without the addition of ferric iron to simulate deep-weathering processes. We report 
here measurements made on the reaction products, and summarize tectonic and geochemical processes that 
might have modified the surface of Mars. 

EXPERIMENTAL PROCEDURES. Powdered samples of pyrrhotite (Fe 7 S 8 ) and olivine (Fa^) were 
reacted in stoppered flasks with sulfuric acid (pH 2) at 60°C, with and without the addition of dissolved 1M 
ferric sulfate, for time periods ranging from 1 to 5 weeks. Reaction products were filtered, washed with 
cold dis till ed water, air-dried and measured by Mossbauer spectroscopy at ambient (295K) and liquid 
helium (4.2K) temperatures to characterize the X-ray amorphous Fe 3+ -bearing phases formed in the 
reactions. Correlative Mossbauer spectra were also obtained for pyrite, marcasite, hematite, goethite, and 
jarosite. The spectra were calibrated against a metallic iron foil standard. 

RESULTS. Examples of 4.2K Mossbauer spectra of pyrrhotite before and after reaction with pH 2 
sulfuric acid are shown in Figure 1. The unreacted pyrrhotite spectrum (Fig. la) is complex and consists of 
3 or 4 superimposed magnetic hyperfine sextets. After reaction with H 2 S0 4 , additional peaks appear at 
approx. 0.1 and 0.7 mm s 1 (Figs lb & lc) attributable to FeS 2 (pyrite or marcasite), which are intensified 
(Fig. lc) when Fe 2 (S0 4 )3 is present as a lixiviant. New peaks also occur near -7.75 and +8.5 mm s" 1 and 
represent another component sextet, the magnetic hyperfine splitting parameter of which (H = 50.6T) 
resembles the value of goethite [5]. The breadth (Fig. lb) or asymmetry (Fig. lc) of the two goethite peaks 
suggest that some hydronium jarosite (H = 49.3T) may also be present in the reaction products. The 
Mossbauer spectra of reacted pyrrhotites, therefore, indicate that sulfuric acid has induced the pyrrhotite -> 
FeS 2 (pyrite or marcasite) -» goethite (+ jarosite) oxidative reactions formulated in Table 1. Nanophase 
goethite is also formed during reactions of olivine with H 2 S0 4 [4], particularly in the presence of dissolved 
Fe 3 * ions, hydrolysis of which is unlikely in acidic (pH 2) solutions. 

DISCUSSION. The experimental results for H 2 S0 4 -degraded pyrrhotites support our proposed model 
for the oxidative weathering of Fe-S sulfide minerals, in which supergene alteration of pyiriiotite in 
groundwater and the formation of secondary FeS 2 phases precedes the deposition of goethite and jarosite in 
gossans above the water table [1-4]. Dissolved Fe 3+ ions clearly promote the formation of supergene pyrite 
or marcasite, and the acidic groundwater causes chemical weathering of ferromagnesian silicates such as 
olivine to nanophase goethite. Similar oxidative weathering reactions involving pyrrhotite and olivine in 
komatiitic basalts are suggested to have occurred on Mars, contributing to the ferric-bearing materials 
present in the regolith there. 

EVOLUTION OF THE MARTIAN SURFACE. Oxidative weathering of sulfides on Mars has been 
largely confined to pyrrhotite, perhaps associated with pentlandite and minor chalcopyrite, which 
predominate as accessory minerals in mafic igneous rocks and in associated ore deposits. Partial melting of 
the martian mantle estimated to contain about 4.5 wt % S [4] produced iron-rich basaltic magma which 
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reached the surface of the planet via massive shield 
volcanoes and as fissure eruptions, transporting high 
concentrations of dissolved S 2 ' and HS to the j M 
surface. Extrusion of turbulent low viscosity lava, 
segregation of immiscible FeS liquids during cooling, 
gravitational settling and fractional crystallization of 
sulfide minerals in the lava flows, produced thinly 
disseminated iron sulfide mineralization over large 
areas of Mars. Massive ore deposits were also 
deposited locally at the base of structural depressions 
or in channels eroded by advancing lava flows [6]. 
Fracturing associated with this volcanism facilitated 
deep-weathering reactions by permeating groundwater 
early in the history of Mars. However, because plate 
tectonic activity appears to have been insignificant on 
Mars, the apparent absence of spreading centers and 94 
subduction zones has minimized acid-buffering of 
aqueous solutions by wall-rock alteration. Therefore, 100 
the acidity of groundwater, now permafrost, has been 
maintained during the chemical evolution of the 
martian surface, aiding the chemical weathering of 
silicates in the basaltic crust when water is present "5 
Furthermore, since tectonic-induced interactions of " 
martian mantle with crust, hydrosphere and 
atmosphere have been minim al, sulfide mineralization 
has not evolved on Mars beyond pyrrhotite- 
pentlandite assemblages associated with mafic 
igneous rocks. Therefore, common terrestrial ore 
deposits such as porphyry copper and molydenum, 
granite-hosted mineralization, and sediment-hosted 
galena-sphalerite assemblages [7] have not formed on, 
Mars. 

References. [1] R.G.Bums, Proc. 17th LPSC, JGR. 92. 
E570 (1987); [2] R.G.Bums, Proc. 18th LPSC, 713 (1988); 

[3] R.G.Bums & D.S.Fisher, LPI Tech. Rept., 88-05 3 4 
(1988); [4] R.G.Bums & D.S.Fisher, LPI Tech Rept., 89-04 ; 
20(1989); [5] E. Murad & J.H.Johnston, in Mossbauer 
Spectroscopy Applied to Inorganic Chemistry (G.J.Long, 
ed.; Plenum)^, 507 (1987); [6] A.K.Baird & B.C.Clark, 
Nature, 311. 18 (1984); [7] C.Meyer, Ann. Rev. Earth 
Planet. Sci., 16. 147 (1988); [8] Research supported by 
NASA grants NSG-7604 and NAGW-1078. 

FIGURE I. Mossbauer spectra measured at 4.2K of pyrrho tite 
(a) before, and (b,c) after, reaction with pH 2 sulfuric acid for 
5 weeks at 60°C, in the absence (b) and presence (c) of 
dissolved 1M ferric sulfate. Peaks attributed to FeS2 are 
identified, as are the positions of the goethite (G) and jarosite 
(J) sextet spectra 
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{1} Fe 7 S 8 + 6Fe 3+ => 4FeS 2 + 9Fe 2+ 

{2} FeS 2 + 14 Fe 3+ + 8 H 2 0 => 15 Fe 2+ + 2 S0 4 2 - + 16 H + 

{3} 2 FeS 2 + 2 H 2 0 + 7 02 (aq) => 2Fe 2+ + 4S0 4 2 - + 4H+ 

{4} 2 FeS0 4 + + FeOH2+ + 6 H 2 0 => (H 3 0)Fe3(S0 4 ) 2 (OH) 6 JJ + 4H+ 
{5} 4 Fe 2+ + 6 H 2 0 + O 2 (aq) => 4FeOOH + 8H+ 

{6} Fe 3+ + 2 H 2 0 => FeOOH + 3 H+ 

{7 } (H 3 0)Fe 3 (S0 4 ) 2 (0H)6 => 3 FeOOH + 2 S0 4 2 ' + 4H+ + H 2 0 
{8} (Mg,Fe) 2 Si0 4 + 4 H + => 2 (Mg 2+ + Fe 2+ ) + SK^ (aq) + 2 H 2 0 
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MOSSBAUER SPECTRA OF OLIVINE-RICH WEATHERED ACHONDRITES: 
II. BRACHINA, CHASSIGNY, ALHA 77005, AND NAKHLA 
Roger G. Burns and Sondra L. Martinez, Department of Earth, Atmospheric and Planetary 
Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139. 

INTRODUCTION. Olivine predominates in several SNC meteorites believed to have originated 
from Mars, including Chassigny (approx. 90 modal % Fa 3 i. 5 ; [1]) and ALHA 77005 (-52% 
Fa28-341 [2]). It is also a significant constituent of Nakhla (15.5% Fa^s; [3]), and is present in 
minor amounts in EETA 79001 lithology A, Lafayette and Govemador Valadares [4]. 
Furthermore, almost 80% of Brachina is olivine (Fa3o) [5,6] (which is related to ALHA 84025; 
[7]), and this meteorite was once grouped with the SNC meteorites. Because olivine is 
particularly reactive in low temperature environments and forms iddingsite under hydrothermal 
conditions and during weathering [8], the presence of such ferric-bearing oxyhydroxide and 
phyllosilicate assemblages may be indicative of oxidative weathering in an aqueous environment on 
Mars. Since the presence of ferric iron in amounts exceeding 1-2% is usually readily detectable by 
Mossbauer spectroscopy [9], and because preliminary measurements of weathered meteorites at 
very low temperatures enabled nanophase hydrated Fe(III) oxide phases to be identified [9,10], 
we have measured the spectra of several olivine-rich achondritic meteorites and report here 295K 
and 4.2K Mossbauer spectral data for SNC and related meteorites; an accompanying abstract [11] 
describes similar measurements on a suite of ureilites found in Antarctica. 

SPECIMENS Studied. Achondrites acquired for the Mossbauer spectral measurements 
included: -100 mg chips from interior portions of Brachina, Chassigny and Nakhla; homogenized 
powders from the Antarctic Meteorite Collection of the "dark" (olivine-rich) and "light" (olivine- 
depleted) fractions of ALHA 77005; and a powdered eucrite (ALHA 80102) showing moderate 
rustiness, which served to locate pyroxene peaks in the low-temperature spectra. Experimental 
details are described in the accompanying abstract [11] 

RESULTS. Assembled in Figure 1 are the 4.2K Mossbauer spectra of the six achondrites. It is 
immediately apparent from Figure 1 that peaks at approx. -7.75 and +8.5 mm s _1 are most 
conspicuous in the Brachina spectrum (Fig. la), yielding a value of 50.3T for the magnetic 
hyperfine field which corresponds to values for goethite. Magnetic ordering of this Fe(III) 
oxyhydroxide phase is not observed in the 295K spectrum, however, indicating that the -7% ferric 
iron present in Brachina occurs as nanophase (<10 nm) goethite. Very weak FeOOH features 
appear in 4.2K Mossbauer spectra of Nakhla (Fig. le) and eucrite ALHA 80102 (Fig. If), but not 
in the spectra of Chassigny (Fig. lb) or the ALHA 77005 "dark" (Fig. lc) and "light" (Fig. Id) 
fractions, even though the 295K spectra resolved small amounts of ferric iron in Chassigny (2.3% 
Fe 3+ ) and ALHA 77005 (1.7% Fe 3 *; [10]) as well as Nahkla (1.6% Fe? + ; [9]). The feme iron in 
Chassigny appears to be accomodated in the olivine structure as it is in ALHA 77005 [12], and 
not in the iddingsite FeOOH phase present in Nakhla and Brachina. 

DISCUSSION. The presence of goethite in Brachina correlates with the "small amount of 
limonite produced by terrestrial weathering" reported when the meteorite was found in South 
Australia [5]. Similar terrestrial weathering could not have occurred in Chassigny and Nakhla 
which were collected as falls, so that their -2.0 Fe 3- *- probably represents preterrestrial oxidation 
during exposure on Mars. We also conclude that the ferric iron measured in olivine-separates [12] 
and bulk samples [9] of ALHA 77005 did not originate in Antarctica but was produced during low 
temperature alteration reactions on the martian surface. 

REFERENCES. [I] RJ.Floran et al., GCA, 42, 1213 (1978); [2] H.Y.McSween Jr et al., Science, 2Q4, 
1201 (1979); [3] T.E.Bunch & A.Miteid, Meteoritics, ID, 303 (1975); [4] H.Y.McSween Jr, Rev. Geophys., 22, 
391 (1985); [5] J.EJohnson et al.,Rec. Sth. Austral. Mus.,11* 309 (1977); [6] C.E.Nehru et al.,Proc. 14th 
LPSC, JGR, £&, B237 (1983); [7] P.H.Warren & G.W.Kallemeyn, Proc. 19th LPSC, 475 (1989); [8] 
R.A.Eggleton, Clays & Clay Min.,21, 1 (1984); [9] T.C.Solberg & R.G.Bums, Proc. 19th LPSC, 313 (1989); 
[10] R.G.Bums, Lunar Planet. Sci., XX, 129 (1989); [11] S.L.Martinez & R.G.Bums, Lunar Planet. Sci., XXL 
this vol. (1990); [12] R.Ostertag et al., EPSL, £1, 162 (1984); [13] Research supported by NASA grants NSG- 
7604 and NAGW-1078. 
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FIGURE 1 . Mossbauer spectra measured at 4.2K of a variety of achondrites. (a) Brachina; (b) 
Chassigny; (c) shergomte ALHA 77005 "dark" (olivine-enriched) and (d) "light" (olivine-depleted) 
fractions, (e) Nakhla; and (f) eucrite ALHA 81002. Vertical dotted lines indicate the locations of 
the outermost peaks of the goethite sextet Mossbauer spectrum. This ferric oxyhydroxide phase is 
present m Brachina, Nakhla and eucrite ALHA 81002. Ferric iron occurring in the olivines, which 
is resolved in the 295K Mossbauer spectra of Chassigny and ALHA 77005, is obscured by peaks 
representing the onset of magnetic ordering of ferrous iron. 
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PHYSIOGRAPHIC CHARACTERISTICS OF MARTIAN DRAINAGE BASINS. CABROL, N.A , 
Lab . Geo . Phys . ( UAO 1 4 1CNRS ) and Lab.Phys.Syst.Solaire.Meudon 92195. -France- 

In this study, a martian drainage basin is defined by the limit with the 
environment of a system of thalwegs. These systems are widly distributed on 
the craterized uplands and display various morphologies and sizes (from 2000 
km 2 to 50 000 km 2 ). A conceptual model, built up with geometrical data from 
the 1: 2.000.000 maps, yields to a classification of drainage basins. This 
classification is based upon the relationship between the perimeter (P) and 
the area (A) of the basin, and is expressed by two coefficients : Kd (drainage 
density) = ZLi/A, where SLi is the sum of the basin branch lengths. Kc 
(drainage basin compacity, Gravellius coefficient) = P/2 tc/A # .28P//A. P is 
the length of a drainage basin limit, which is defined in this method by a 
line located 10 km upstream the headwater system and including the branch 
systems until the outlet. This 10 km value considers the limited Viking 
resolution, and the limit of the hydrogeological system of underground water 
reservoirs. To test this model, 15 well preserved drainage basins are selected 
between -47” and +19”Lat. in order to cover a representative distribution. The 
model is calibrated by 8 terrestrial basins, located on 3 different geological 
units (Bassin Aquitain, Bassin de la Durance and Bassin de la Seine -France-), 
respectively in 3 rainfall regimes. These coefficients Kc and Kd characterize 
the basin physiographies. Within the model, a discriminant parameter Kd/Kc = D 
(drainage intensity) is established. Terrestrial Kc values vary from 1 to 1.5 
and Kd from .02 (karsts : scarcely branched areas) to 300 (badlands, high 
drainage intensity areas). So, terrestrial D is included between 200 (high 
drainage intensity) and .01 (low drainage intensity) with medium compacity 
value of 1.5. The relationship between Kd and Kc allows to classify martian 
and terrestrial basins (Fig.l) which are represented on the graph by their 
parameters Kd and Kc. Three concentrations appear, one for the terrestrial 
examples, and two for the martian basins. The gravity center for the 
terrestrial examples is Kd= .05 and Kc = 1.3, so D = .04. For the martian lower 
concentration I, the gravity center coordinates are : Kd=.l and Kc=1.6, so 
D=.06. For the martian concentration II, the gravity center coordinates are : 
Kd= . 05 and Kc=3, so D=.01.The D martian values for the 15 tested basins are 
closer to the representative value of terrestrial low drainage basins 
intensity located on geological units of karstic types. The martian drainage 
basins of (5 + /-2) .10® km 2 correspond to a high compacity value and a low 
thalweg distribution. Martian drainage basins of (30 + /-10) . 10® km 2 correspond 
to a low compacity value and a low thalweg distribution. 

A first conclusion is that Kd value varies only between .03 to .15 (Earth .02 
to 300). The result is calibrated with 30 other drainage basins on Mars. 
Because Kd is related to the geological type of the basin, this result 
indicates a probable undiversified substratum on a wide area (-50”, +10”Lat.) 
in which local singularities might lead to groundwater occurrences. A second 
conclusion is : for martian drainage basins, Kc appears as a clear 
discriminant coefficient. 

These conclusions are then correlated with a mapping of headwater branch 
densities (order 1) on the entire martian surface. The results vary from 100 
to 1 per 10® km 2 , implying a wide interval of variation. These headwater 
systems are found mainly in the Margaritifer Sinus, Sinus Sabeus, and Iapygia 
regions. Their dispersion is unconsistant with an organized meteorological 
feeding origin and with true surface drainage basins. The .03 mean Kd value 
and the 1.6 to 3 Kc values (corresponding to elongated basins with isolated 
headwater systems) lead to support the seepage hypothesis. 
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The third conclusion is : the limits of martian drainage basins are not 
divides. The basins are not joined, they are widely spreaded as relatively 
small units on wide areas. The methodological hypothesis considering these 
basins as isolated systems is consistant with the observation. 

REFERENCES : Abrahams, A.D, 1984 Channels networks : a geomorphological 
perspective. Water. Res. Research, 20 , pl61-188. Baker, V.R, and Partridge, J.B, 
1986 Small martian valleys : Pristine and degraded morphology. J.G.R 91, 
p3561-3572. Battistini, R, 1985, in Klinger et al (eds), Hydrolithosphere and 
problem of subsurface ice in equatorial zone of Mars, Ices in the Solar 
System, p607-617. Cabrol, N.A, 1988, Compte-rendu des Journees de 
Planetologie. (CNRS/INSU). Cabrol, N.A, 1988 Morphological variations and 
evolution of channels on Mars. LPI , XX, pl36-137. Carr, M.H, 1986 Mars a water 
rich Planet? Icarus 68, pl87-216. 
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ANALYSIS OF MARINER 6 AND 7 SPECTRA FOR WEAK ABSORPTION FEATURES FROM 2 
TO 6 pm, Wendy M. Calvin and Trude V. V. King, U. S. Geological Survey, 
Denver, Colorado 80225 


The current study of the Mariner 6 and 7 spectra was undertaken to iden- 
tify previously unrecognized absorption features and to confirm the results 
of recent ground-based observation of Mars by Clark et al . ( 1989 ), Blaney and 
McCord (1989), and Pollack et al. (1989). The Mariner data have the advan- 
tage that they are not affected by terrestrial atmospheric absorption, have 
higher spatial resolution than terrestrial observations, and cover the short 
wavelength infrared region with reasonably high spectral resolution. How- 
ever, the Mariner data are subject to a variety of other uncertainties which 
we have attempted to correct for, in as far as is possible twenty years after 
the original flyby. 

The Mariner 6 and 7 infrared spectrometer data set was obtained from 
Martin (1985), who provided corrections to the current latitude and longitude 
system and an initial wavelength calibration. The wavelength region from 2 to 
6 pm is covered by two halves of a circular variable filter wheel. The first 
wavelength segment is 1.9 to 3.7pm and the second is 3.0 to 6.0pm. In 
between each half of the filter wheel is a clear aperture which produces a 
fiduciary spike ostensibly for wavelength calibration. The use of the fidu- 
ciary spikes for this purpose was found to yield inaccurate wavelengths, as 
did corrections using atmospheric bands at one side of a wavelength segment. 
Apparently the rotation rate of the filter wheel varied nonlinearly and so 
adjustments at one side of the wavelength segment often will not accurately 
predict wavelengths on the other side of the wavelength segment. Further 
wavelength calibrations were performed using atmospheric CO^ bands, where 
available, and fiduciary spikes. In the 3 to 4-pm wavelength region where no 
strong atmospheric bands occur a feature caused by an absorption in the 
filter at 3.4pm (Herr et al. , 1972, and Pimentel et al . 1974) was used for 
calibration. 

Other corrections made to the data set include removal of a space back- 
ground signature, correction for the instrument response as a function of 
wavelength, removal of the thermal flux to the reflected component and 
finally correction of the reflected solar response to reflectance. These 
additional corrections were only made for the 3 to 6- fj.ru wavelength region 
because at shorter wavelengths the space background did not vary signifi- 
cantly from zero and calibration spectra are not yet available to derive 
accurate instrument response functions. For the longer wavelength segment 
space background responses were determined from 6 spectra recorded after the 
flyby. The average of these was fairly noisy and it was necessary to smooth 
the response function using a fast Fourier transform. Instrument response 
functions, determined from black-body calibration^ spectra included in the 
flight data set, were also smoothed to eliminate the introduction of noise. 
In addition, the solar flux was calculated and adjusted for Mars' distance 
from the sun, and the thermal contribution was calculated using temperatures 
determined by the original instrument team (Pimentel et al . , 1974). 

In the wavelength region from 1.9 to 3.6pm, several absorption features 
that cannot be attributed to atmospheric absorptions are observed in the pro- 
cessed Mariner 6 and 7 data. The most dominant feature is the broad 3pm 
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absorption associated with hydrated surface minerals which was discussed by 
Pimentel et al . (1974). An absorption feature is also observed in the 2.3 to 
2.4-jzm wavelength region which is in agreement with the Clark et al . (1989) 
observations. However, the resolution of the instrument is insufficient to 
verify whether the 2.3-2.4^m feature is atmospheric or mineralogic in origin. 
Between approximately 2.45 and 2.62pm there are two broad absorption bands. 
The shorter of the two bands is centered on 2.48^m and seems to be correlated 
to a band observed by Clark et al. (1989) at 2.45/im. The 2.48//m absorption 
feature in the Mariner data could be shifted in wavelength because the 
wavelength calibration is fixed at isolated points, the closest of which is 
at 2.7yum. The second band is broader (2.5- 2.65^m) and can possibly be 
explained by discrepancies between the amount of atmospheric water modeled 
and observed or the presence of hydrated surface minerals. Removal of the 
atmospheric signature in the 2.7 to 2.8-pm region indicated a residual 
absorption feature at 2.76/im which can be attributed to the presence of OH- 
bearing mineral phases. This feature can be seen weakly in some of the origi- 
nal Mariner 7 spectra. Attempts to confirm it with Mariner 6 were unsuccess- 
ful as the signal to noise ratio of Mariner 6 data is much worse than that of 
Mariner 7. Uncertainties associated with removing atmospheric absorptions 
have precluded a definitive identification of an OH-absorption feature. 

In the 3 to 6-pm region the strongest absorption features result from 
atmospheric absorptions. However, a broad 3-/im absorption associated with 
hydrated minerals is present. The 3.8-jum feature reported by Blaney and 
McCord (1989) is confirmed. Over Hellas the 3 . 8^m feature is nearly twice as 
deep as is predicted for a four airmass atmosphere, indicating a mineralogi- 
cal as well as atmospheric origin. In addition, a weak feature at 5.35/im, on 
the edge of the atmospheric CO^ band at 5.2/im, is observed in both Mariner 6 
and Mariner 7 spectra of the Meridiani Sinus region. These absorption 
features can result from the presence of OH-bearing mineral species and pos- 
sibly, based on preliminary laboratory spectra, from S-bearing mineral phases 
or a combination of the two mineral species. 

This study confirms the presence of several weak absorption features 
observed by terrestrial observations and identifies previously unrecognized 
absorption features in the wavelength regions from 2.5 - 2.65 and 5.35^m in 
Mariner 6 and 7 data. These features appear to be mineralogic rather than 
atmospheric in origin and are believed to result from the presence of OH- or 
S-bearing minerals on the Martian surface. 
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THE EFFECTS OF FLOODS, VOLCANISM AND POLAR PROCESSES ON THE D/H RATIO 
IN THE MARTIAN ATMOSPHERE. M. H. Carr, U. S. Geological Survey, Menlo Park, CA 
94025. 


Water in the martian atmosphere is 5.1 times more enriched in deuterium 
than terrestrial water (Owen et al. , 1988; Bjoraker et al., 1989). The 
enrichment has been previously attributed to either massive loss of water early 
in the planet's history (Owen et al., 1988) or to the presence of only a very 
small reservoir of water that has exchanged with the atmosphere over geologic 
time (Yung et al. 1988). However, early massive loss of nearly all the planet’s 
water appears inconsistent with abundant evidence for the action of water 
throughout Mars history, and exchange of only small amounts of water with the 
atmosphere is contrary to interpretation of channels as water worn, and 
interpretation of the poles as a sink for water. Here we examine what effect 
geologic processes might have on D/H in the atmosphere, and assess whether the 
geologic model for the evolution of the Mars surface that has been developed over 
the last two decades needs major revision. 

Outflow channels are generally believed to have been formed by large 
floods. A large flood would immediately inject a substantial amount of water 
into the atmosphere, then in the long term water would sublime into the 
atmosphere from the ice rich deposits that must have formed where the water 
pooled after the flood was over. To estimate the amount of water injected into 
the atmosphere, a flood was divided into three stages. Assuming present 
conditions on Mars at the time of a flood, and using previously derived 
techniques for determining sublimation rates (Toon et al. , 1980), it was 
estimated that roughly 10 17 g of water would evaporate into the atmosphere in 
the active stage, while the flood was in progress. This water probably would 
have frozen out rapidly at the poles to form an ice layer a few centimeters 
thick. In the second stage, the lake that formed at the end of the flood channel 
freezes, and the thermal anomaly created by the presence of the lake dissipates. 
It was esti ma ted that in a typical flood the terminal lake would take about one 
year to freeze solid and the thermal anomaly would take about 6 years to 
dissipate during which time about another 10^ gm of water would evaporate into 
the atmosphere and be precipitated out at the poles. The long term fate of the 
terminal ice deposits would depend on their latitude and thickness of any 
superimposed dust. Most of the floods in Elysium and Hellas, and around Chryse 
Planitia appear to have pooled at high latitudes (McGill, 1985; Lucchitta et al. , 
1986). Dust would affect the stability of ice deposits by reducing daily 
temperature excursions and inhibiting movement of water vapor from the buried 
ice to the surface. These effects were modelled using the techniques developed 
by Clifford and Hillel (1983). For an ice deposit at 50° latitude, covered with 
a few tens of centimeters of dust, sublimation rates are estimated to be 
extremely low, in the range of 10 -7 - 10” g cm yr . Sublimation rates at 
the equator are over 100 times higher. Thus, the immediate effect of a flood on 
the atmospheric D/H would be to reset the atmospheric D/H to the value for 
groundwater. The long-term effect would depend on the latitude of the terminal 
ice deposit. Sublimation rates are so low for high latitude deposits that the 
D/H in the atmosphere could become enriched with respect to groundwater in 
response to upper atmosphere losses. If the terminal ice deposit were at low 
latitude, the atmosphere would be maintained at the D/H value for groundwater 
until the deposit had fully sublimed. 
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Volcanism is likely to have episodically introduced water into the 
atmosphere throughout the history of the planet. Greeley (1987) estimates that 
25 to 41 m of water averaged over the whole planet have been introduced onto the 
surface by volcanism during the last 3.5 Ga. Judging from the size of lava 
flows, individual eruptions could almost instantaneously introduce 10 13 to 10 17 
gm of water into the atmosphere. The larger eruptions would inject enough water 
into the atmosphere to reset its D/H ratio to that of juvenile water. Much of 
the erupted water would have rapidly frozen out at the poles without having its 
D/H ratio changed. 

The effect of polar processes on D/H in the atmosphere depends on how 
actively water at the poles is exchanging with the atmosphere. If the Yung et 
al. (1988) estimates of loss of hydrogen from the upper atmosphere are correct, 
and water at the poles has actively interchanged with the atmosphere, then the 
poles can contain no more than 0.2 m of water averaged over the whole planet. 
But this conclusion is inconsistent with evidence of introduction of water onto 
the surface by floods and volcanism, and with the current perception of the 
polar layered terrains as ice rich. Jakosky (in press) suggested that Yung et 
al (1988) had underestimated hydrogen loss rates. He suggested that loss rates 
at high obliquity could be 100 times present rates, and the time averaged value 
20 times the present rate. If so, then the reservoir of water at the surface 
would be correspondingly increased. An alternative explanation of the 
inconsistency is that exchange of water between the poles and the atmosphere is 
minimal, and restricted to the upper few centimeters of the polar deposits. 
Thus, new ice added to the polar deposits as a result of floods or volcanism, 
retains its initial D/H ratio except for the upper surface. In this case, the 
D/H in the atmosphere evolves almost independently of the poles, and reveals 
little about the inventory of water at the surface. Such a scenario is 
consistent with the crater ages of Plaut et al. (1989), which suggest that the 
polar layered terrains are older than formerly thought. 

We can conclude that the D/H enrichment of the atmosphere is not 
necessarily that of the bulk of the near surface water. The evolution of D/H 
in the atmosphere is likely to have been episodic. Floods and eruptions would 
have periodically reset the atmospheric D/H to the groundwater and juvenile 
values respectively. Between these events the atmosphere would have become 
enriched in deuterium at a rate dependent on the rate of interchange of water 
between the atmosphere and other surface reservoirs such as the poles. 

REFERENCES. Bjoraker, G. 1., Mumma, M. J., and Larson, H. P. , Proc. 4th Int. , 
Conf. on Mars, Tucson, 69-70 (1989). Clifford, S. M. , and Hillel, D., J. Geophys. 
Res., 88, 2456-2474 (1983). Greeley, R. , Science, 236 , 1653-1654 (1987). 
Jakosky, B. M. , J. Geophys. Res. (in press). Lucchitta, B. K. , LPSC XVII, 498- 
499 (1986). Me Gill, G. E., LPSC XVI, 534-535 (1985). Owen, T. , Maillard, J. 
P., deBergh, C., and Lutz, B. L. , Science, 240 , 1767-170 (1988). Plaut, J. J. , 
Kahn, R. , Guiness, E. A., and Arvidson, R. E., Icarus, 76 , 357-377 (1988). Toon, 
0. B., Pollack J. B, Weird, W. , Bums, J. A., and Bilski, K. , Icarus, 44 , 552- 
607 (1980). Yung, Y. L. , Wen, J. , Pinto, J. P. , Allen, M. , Pierce, K. K., and 
Paulsen, S. , Icarus, 76, 146-159 (1988). 
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CRATER MORPHOLOGY VARIATIONS IN THE ELYSIUM REGION: 
IMPLICATIONS FOR ICE DISTRIBUTION ON MARS 

Julie A. Cave, U.L.O. Planetary Image Centre, 33-35 Daws Lane, Mill Hill, 

London NW7 4SD. 

INTRODUCTION: Variations in crater morphology are attributed to differences in target strength, 
depth of excavation, presence of near-surface volatiles, and environmental conditions (1,2, 3, 4). Rampart 
craters and other features which may be a result of the presence of sub-surface ice have been observed in 
the Elysium Volcanic province (5). The aim of this study is to examine the distribution of various crater 
characteristics and to evaluate the usefulness of crater morphology as an indicator of sub-surface ice. 

METHODOLOGY: Preliminary mapping of the region indicated a great complexity of crater mor- 
phologies. The study area (155-245°, 15°S-47.5°N) covers a wide range of latitude, altitude, surface type 
and age and there are therefore several factors that could have contributed to the observed morphological 
variations. Using mosaics of Viking frames (at a scale of 1:1250000) the co-ordinates of each crater over 
1.875 km in diameter were digitised. A comprehensive classification scheme was devised that recorded the 
characteristics of each component of the crater allowing as much detail as possible to be used in the anal- 
ysis. The characteristics of the craters’ morphology recorded include: profile, central features, rim nature 
(continuous, buried etc.), rim condition (fresh, terraced, degraded efc.), ejecta type and surface texture 
(for each component of multiple blanket ejecta), geological unit according to references (6) and (7), and 
local target nature ( e.g . lava flow, crater ejecta, channel floor). A confidence factor was also included since 
the precision of interpretation varied with resolution, although wherever possible the descriptive details 
were confirmed by checking the highest resolution images available. These details were incorporated into 
a homogeneous database which can be searched systematically to highlight the dominant factors. Initial 
analysis has concentrated on average or dominant crater characteristics as a function of latitude, altitude, 
geological unit, crater diameter and local target nature. The database is nearing completion, and initial 
observations presented here come from a preliminary analysis of 3060 craters in the Northern Hemisphere. 
OBSERVATIONS: 

(a) Approximately one third of the selected craters have discernable ejecta blankets. Plots of maximum 
ejecta diameter versus crater diameter (hereafter abbreviated to Evs$) provide an indication of the mobility 
of the ejecta which is presumed to reflect the degree of fluidisation. The majority of the ejecta documented 
showed some degree of fluidisation (even for the smallest craters), and very few radial ejecta blankets were 
observed (this may be due in part to the resolution, since radial deposits are harder to detect). If the 
fluidisation of the ejecta is due to the incorporation of subsurface volatiles, then variations in the gradients 
of Evs$ reflect changes in the state or concentration of volatiles rather than a transition from dr> to 
‘wet’ materials being excavated (ignoring possible atmospheric effects). 

(b) For most geological units there is a well-defined relationship between ejecta diameter and crater 
diameter. Many units (although the statistics of some are limited) seem to exhibit significantly higher 
gradients at larger crater diameters. When all the data are combined the resulting graph is a smooth, 
rising curve, but the data for several units taken individually exhibit gradient discontinuities. The position 
of such discontinuities, and the gradients appear to depend on the geological unit. The Elysium lava 
plains, unit AEL1 have a lower gradient of 2.1±0.2, upper gradient of 3.6±3.6, with an assumed break- 
point occurring at crater diameter 7km. The lower gradient lies between the gradient for similar-sized 
craters on Elysium Mons itself (AEL2, 2.0±0.1) and that of the channeled deposits AEL3 which originate 
from the N.W. flanks of the volcano (gradient = 2.9±0.2). This might suggest a localised transfer of 
volatiles from AEL2 to AEL3 driven by the increased thermal gradient of Elysium Mons. Further evidence 
of volatile enrichment of units AEL3/4 is provided by the marked concentration of double-ejecta craters 
in this location (see figure 1). 

(c) An initial investigation of Evs$ has revealed no obvious trend with latitude, except that the 
previously mentioned relationship breaks down at latitudes greater than about 40° N. Plots of Evs4 in 
this region indicate a wide scatter of the data, which may be due to the high concentration of pedestal 
craters here. 

(d) The precise target nature does appear to have influenced the spread of the ejecta: a linear Evs$ 
relationship occurs for craters that have impacted other crater rims (despite varied geological locations) 
and a possible linear trend has been observed for craters on volcanic constructs. 
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(e) The location of central features is as shown in Figure 1. The peaked craters are widely distributed 
on many units, while the pit craters are less frequent and have a more limited distribution. The plot of 
Evs$ for these craters indicate similar ejecta fluidities, possibly indicating similar formation conditions. 
The scarcity of pits may be due in part to the infilling or flooding of many of the larger craters in the 
sample set. 

(f) These observations are currently being checked for statistical accuracy. Altitude, radial distance 
from Elysium Mons, and estimates of the depth of excavation will also be taken into account. The 
morphological characteristics of the ejecta will be examined simultaneously with the dimensional data. 

The database is a valuable aid in the interpretation of the crater morphology in the area, and will 
also provide information on the ages and resurfacing histories of the units. This work forms part of a 
detailed investigation of Elysium in which all possible indications of the presence of ice (craters, channels, 
chaotic and knobby terrain, volcano eruptive styles) are being assessed. Comprehensive studies of the age, 
nature and distribution of each ice-related landform will enable further constraints to be placed on the 
importance, timing and distribution of ice in Elysium. In a future project, the Hellas volcanic region will 
be investigated using the same technique, allowing comparisons to be made between two volcanic provinces 
of different ages and geological settings. 

REFERENCES: (1) Head J.W., Proc. L.P.S.C. VII 2913-2927, 1976. (2) Carr M.H., Crumpler L.S., 
Cutts J.A., Greeley R., Guest J.E. and Masursky H., J. Geophys. Res. 82 4055-4065, 1977. (3) Mouginis- 
Mark P.J., J. Geophys. Res. 84 8011-8022, 1979. (4) Bridges N.T. and Barlow N.S., L.P.S.C XX 105-106, 
1989. (5) Mouginis-Mark P.J., Wilson L., Head J.W., Brown S.H., Lynn Hall J. and Sullivan K.D., Earth, 
Moon, and Planets 30 149-173, 1984. (6) Scott D.H. and Tanaka K.L., Map I-1802-A, U.S.G.S., 1986. (7) 
Greeley R. and Guest J.E., Map I-1802-B, U.S.G.S., 1987. 
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CROSSCUTTING, PERIODICALLY SPACED WRINKLE RIDGES OF HESPERIA PLANUM; 
D. John Chadwick, Dept, of Geology, Northern Arizona University, Flagstaff, A Z 86001 ; 
Thomas R. Watters and Michael J. Tuttle, Center for Earth and Planetary Studies, National 
Air and Space Museum, Smithsonian Institution, Washington, D.C. 20560 


Wrinkle ridges are the long, sinuous anticlinal features which occur on the volcanic 
plains units of the terrestrial planets and the Moon. Most of these ridges form single, 
parallel trending groups or concentric patterns within circular basins. The wrinkle ridges 
of Hesperia Planum on Mars, however, show a much more complex spatial relationship. 
Two distinct sets of ridges have formed here, crosscutting one another at nearly 
orthogonal angles, implying a more complex structural history for the region. Termed 
reticulate ridge patterns (1), these crosscutting ridges are found elsewhere on Mars, but 
at Hesperia Planum they dominate a very large area to the east and southeast of Tyrrhena 
Patera. 

In this study, the two directional components of the reticulate pattern were 
separated and ridge spacing was determined using a series of sampling traverses spaced 
about 20 km apart and oriented perpenducular to the mean ridge trend of both sets. In 
addition, the two sets of ridges maintain a constant trend over a limited area. Therefore, 
a large area of Hesperia Planum containing the reticulate pattern was divided into four 
domains based on the mean trend of the ridges. In each of the domains, the ridge 
spacings of the two sets were found to be very similar and passed the statistical t-test for 
equivalency. This test compares two sets of sampled numerical data to determine if they 
appear to be derived from the same population. The mean spacing for the two ridge 
sets in each domain are as follows: 38.23 km and 36.13 km for domain one, 32.02 km 
and 33.75 km for domain two, 36.10 km and 35.58 km for domain three, and 23.80 km 
and 24.26 km for domain four. The equivalent spacing of both sets of ridges can be 
explained by two superimposed episodes of buckling at a critical wavelength, governed 
by the rheological properties and thickness of the plains material in which the ridges 
formed (2). 

The origin(s) of the compressive stresses which formed the reticulate pattern is not 
readily apparent. The compressional stresses that resulted in the formation of the wrinkle 
ridges of the Tharsis Plateau may have been in part due to isostatic uplift (3, 4,5, 6). In the 
absence of a “Tharsis-like" uplift or load in Hesperia Planum, other mechanisms for the 
generation of compressional stresses must be found. Compression may have resulted 
from subsidence due to loading of the volcanic plains material (1), but this mechanism 
may produce only a single ridge trend related to the shape of the basin. A superimposed 
regional tectonic event may account for the second ridge set in the reticulate pattern, with 
local influences accounting for the variable trends of the ridges. 
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CROSSCUTTING WRINKLE RIDGES: D.J. Chadwick et al . 



Figure 1. Crosscutting 
wrinkle ridges of Hesperia 
Planum; a portion of domain 
three. 


1 - 6 $^' 



100 km 


Figure 2. The four domains of 
eastern and southeastern Hesperia 
Planum. 
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GEOLOGIC MAPPING OF LOWER MANGALA VALLES, MARS: EVIDENCE OF FLOODING, 
SAPPING, DEBRIS FLOW, AND VOLCANISM 1 ; Mary G. Chapman and Kenneth L. 

Tanaka, U.S. Geological Survey, 2255 North Gemini Drive, Flagstaff, AZ 
86001. 

Two periods of catastrophic flooding are indicated by detailed 
geologic mapping of the Mangala channel, based on two 1: 500,000-scale 
photomosaics (MTM -05152 and -10151), computer enhancement of Viking 
images, and photoclinometric channel profiles. The ages of these two 
episodes are Late Hesperian and Amazonian-Hesperian. Lobate lava flows 
were erupted during channeling episodes; at least two stratigraphically 
distinct, resistant lava flows filled the early channels; later flooding 
eroded the surrounding, less resistant material to leave the lava-capped 
channels higher than surrounding terrain (i.e. inverted topography). 

Impact craters that were breached and filled by the earliest phase of 
flooding formed lakes; some of these ponded areas were covered by lava, 
others were refilled by the later period of flooding. This flooding period 
was followed by (1) spring sapping that formed narrow, theater-headed 
channels (2) debris-flow deposition resulting from chaotic collapse of 
lava-capped material within breached craters. The direction of drainage 
was north toward Amazonis Planitia, where the flow materials were partly 
covered by Middle Amazonian lobate lava flows. 


‘Previously published in GSA Abstracts with Programs, 1989, v. 21, no. 6, 
p. 4108. 
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ESA PLANS FOR THE SCIENTIFIC EXPLORATION OF MARS 

Chicarro A.F., Space Science Department, ESA/ESTEC, 2200 AG Noordwijk, The Netherlands. 


The exploration of Mars by unmanned spacecraft during the last quarter of a century shows that, 
although geologically less evolved, Mars is far more Earth-like than any other planet. Its internal 
evolution, resulting in intense tectonic and volcanic activity, and the exogenic processing of its surface, 
have extended over several billion years. Aside from the Earth, Mars is the only other planet with 
surface-temperature conditions in the range of stability of complex organic compounds, and therefore, 
it is an obvious target in the search for present or extinct lifeforms. 

Specific goals: The future international exploration of Mars must address major scientific questions, 
such as: What is the present internal structure, composition and activity of the planet? What are the 
relative and absolute ages of Mars geologic units? What is the chronology of extensive and compressive 
tectonic phases? Do SNC meteorites unequivocally come from Mars? What are the most common 
minerals in volcanic, plutonic and sedimentary rocks? How different were the atmospheric composition, 
pressure, temperature and dynamics in the past? What is the present surface and subsurface water 
distribution ? How closely did early Mars resemble the Earth when life appeared on our planet, about 
3.8 billion years ago? Are there sediments on Mars containing organic compounds or primitive fossils? 

Strategy: ESA has identified three possible areas of European participation in the future international 
exploration of Mars. These options, two of them dealing with science from the surface, the other 
addressing science from orbit, are as follows: 

• Network science: A network of small surface stations, including a combination of 2 hard landers 
(surface penetrators) and a cluster of 3 semi-hard landers (mini-probes), thus defining a global/local 
seismological network (Figure 1) designed for an operational lifetime of one Martian year. The scientific 
objectives (Table 1) would especially concentrate on the internal structure, mineral and chemical 
composition, surface meteorology, and entry and descent science measurements. The global/regional 
network concept could be part of a precursor mission to the Mars Rover Sample Return and Manned 
Exploration missions. 

• Rover science: A sophisticated sample-acquisition subsystem, including intelligent robotic arm(s) to 
be placed on-board a rover from another space agency, and capable of stereo- vision, handling, sampling 
and analyzing rocks and soils. It would conduct detailed in-situ chemical and mineralogical analyses 
of surface samples, subsurface structural studies, biological experiments, and atmospheric studies, at 
a number of surface locations within a limited area of mobility. This sampling subsvstem could be a 
contribution to a Mars Rover Sample- Return mission. 

• Orbiter science: A facility on-board an orbiting spacecraft able to: (i) acquire very high-resolution 
imagery, (ii) provide radar altimetry and subsurface sounding, or (iii) contribute to the identification 
of the mineralogical and chemical composition of geologic units. This orbital facility would include 
a very sophisticated imager or a microwave radar instrument, and would have scientific potential for 
any surveillance mission of Mars. 

International cooperation: Any of the three selected areas for possible European participation 
would represent a major and independent contribution to a future international mission to Mars, 
complementing the already scheduled US Mars Observer and the USSR’s Mars-94 missions. However ^ 
in-depth technical and scientific assessment studies, with emphasis on a Network mission, are currently 
under way in ESA. Following the cooperative path initiated by Cassini, a participation in international 
planetary missions stands as a sensible and viable approach for Europe, in order to take its proper 
role in the exploration of the Solar System during the next decade and beyond. 

Reference. Chicarro A.F. et al., Mission to Mars: Report of the Mars Exploration Study Team, ESA 
SP-1117, European Space Agency, 138 pp., 1989. 
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Figure 1: Location of proposed landing sites in the Tharsis region. 


Site ' Device 

1 Cluster of 3 mini- 
probes 


Total of: 

I global seismometer, 
3 local seismometers 
(1 axis), 

3 meteorol. packages, 
1 camera, 
chemical analysis 
(DEG, TPA, PM, 
GC)' 


Independently targetted 
penetrator 

For each penetrator: 
global seismometer 
(3 axis). 

Independently targetted camera, 
penetrator meteorol. package, 

chemical analysis 
(DEG, TPA, XRF, 


250 kg including 
aeroshell and 
- 5 kg payload on 
each probe 


100 kg including 
- 8 kg payload 


100 kg including 
- 8 kg payload 


Geochemistry 

Analyse ridged plains 
material and subsur- 
face volatile content 


Scientific interest 
Meteorology 

Northern hemisphere 
site, near Viking 
Lander- 1 in order to 
compare results 


GRS, GC, PM, CS) 


Analyse most wide- 
spread older unit: 
complex history of 
weathering and im- 
pacting 

Analyse younger 
volcanic unit: type of 
volcanic activity. 

All three sites have 
simple geology: 
unmistakable units 


Southern hemisphere 
site: study formation 
and evolution of 
global dust storms 


Near the equator and 
higher-altitude site: 
study formation of 
clouds in addition to 
weather patterns and 
atmospheric circulation 


Seismology 

The three sites deter- 
mine a triangle (of 
side ~ 3500 km) 
which englobes most 
of the Tharsis uplift 
and its radial fault 
system: 

— likely tectonic ac- 
tivity 

— good definition of 
seismic epicentres. 
Global network: to 
determine Mars’ inter- 
nal structures 

Local network: to 
study the lithosphere’s 
upper layers. In addi- 
tion, meteoritic im- 
pacts do trigger 
seismic waves: to help 
determine the planet’s 
internal structure 


Table 1: Scientific interest of a network of small stations. 
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DERIVATION OF MARS ATMOSPHERIC DUST PROPERTIES FROM 
RADIATIVE TRANSFER ANALYSIS OF VIKING IRTM EMISSION PHASE 
FUNCTION SEQUENCES; R.T. Clancy and S.W. Lee, Laboratory for Atmospheric and 
Space Physics, University of Colorado, Boulder, CO 80309 


During the span of the Viking Orbiter missions, several hundred emission-phase-function 
(EPF) sequences were obtained, in which the IRTM instrument observed the same area of surface 
as the spacecraft moved overhead. The IRTM data set as a whole is well calibrated, having been 
con-ected for inter-spacecraft, inter-detector, and temporal calibration variations; a conservative 
estimate for the absolute uncertainty inherent in this data set is 1-2% [1]. The EPF visual 
brightness observations (passband of 0.3-3.0 pm , effective wavelength of 0.57 pm ) present 
veiy accurately calibrated albedos of a given region and the atmosphere above it versus emission 
(and to a lesser extent, incidence) angle. Although these observations were designed for the 
purposes of surface photometric studies, they are quite suitable for the derivation of atmospheric 
dust opacities. 

We have developed a radiative transfer model of the Mars atmosphere and surface based 
upon a discrete-ordinates radiative transfer code [2]. Input parameters include the atmospheric 
dust opacity, the single scattering albedo (GJn) and particle phase function of the atmospheric dust, 
and the surface bidirectional reflectance. A 16-stream approximation is adopted to adequately 
represent the model dust and surface phase functions. The surface photometric function is 
estimated from EPF sequences at low dust opacities and emission angles, where the effects of 
scattering by atmospheric dust are minimized. The exact behavior of the bidirectional reflectance 
at large emission angles and/or high dust loading is not critical since the contribution of 
atmospheric dust scattering dominates the observed reflectance for these regions. 

On the basis of an analysis of ~20 separate EPF sequences spanning low and high latitudes, 
low and high dust opacities, and a wide range of L s ; we find evidence for three distinct aerosol 
populations. In all three cases it is possible to model the dust scattering with a single dust phase 
function labelled "IRTM" in figure 1 . An empirical phase function representing scattering by 
non-spherical particles was derived by Pollack et al. [3] to fit the Viking lander observations of 
dust scattering, and is presented by the solid line in figure 1. The agreement between the "IRTM" 
and Pollack phase functions is quite good; however the asymmetry factor, g, is 0.55-0.56 rather 
than the value of 0.79 given by Pollack et al. [3]. A Henyey-Greenstein phase function with 
g=0.79 is also included in figure 1 to indicate the very large large differences in particle back- 
scattering implied by g=0.55 versus g=0.79. 

We find three separate values of the single scattering albedo, GJ 0 , corresponding to the three 
aerosol populations. At low-to-mid latitudes for moderate dust opacities (t = 0.2- 1.0), G5 0 =0-88 
fits die EPF observations. Figures 2-4 present a subset of model-data comparisons for this case. 
At high latitudes (70°-80°N) we find B3o=1.0, indicating that ice clouds are the likely source of 
particulate scattering. Figure 5b presents a model-data comparison for this case and a best fit 
attempt with GJq= 0.88 (figure 5a). The third aerosol population corresponds to very high dust 
loading (x > 2), in which it is necessary to increase GJq to 0.9-0.92 to fit the data with the "IRTM" 
dust phase function. Figures 6a and 6b indicate model-data comparisons over Hellas basin in 
which we attempt to fit the EPF data with B5 0 =0.88 and G5q=0.92, respectively. 

The primary differences between the EPF IRTM and Viking lander derived dust properties 
are larger GJ 0 from the EPF analysis (GJo=0.88 versus 0.86) and an increase in G$o during the peak 
of the dust storm as found in the EPF analysis (the Viking lander data indicate no such change, J. 
Pollack, personal communication). It is possible that the distinction of looking from below with 
lander data and looking from above with the IRTM data may account for these differences. Finer 
dust (higher GJq and more backscattering) may be mixed to higher altitudes and may increase 
during the peak of the dust storms (J. Bergstrahl, personal communication). The downward- 
looking geometry of the EPF sequences would be more sensitive to a high-altitude, fine 
component. 

This research was supported under NASA Planetary Geology grant NAGW 1378. 
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COLD POLAR DESERTS: WEATHERING RATES, MECHANISMS, AND IMPLICATIONS FOR REMOTE 
SENSING OF MARS; Edward A. Cloutis, Department of Geology, University of Alberta, Edmonton. Alberta T6E S6 

Introduction- Polar and alpine desert regions (cold deserts) offer the best terrestrial field analogues to^ weathering 
conditions on Mars [1]. A suite of P rock samples from the Northwest Territories in Canada has been s^tral y ch^actemed 
in order to better understand the weathering processes prevalent in cold deserts and to dadlri onto 
sensing of Mars. Since these samples have presumably been subaenally exposed only since Lf ' * 

are useful for exami ning how relatively short-term physical and chemical weathering affect spectral properties. 

USC Thesamph^were collected from the Yellowknife Bay area (basalt) and Contwoyto Lake-Pomt Lake ; region m th 
Slave Province. The samples include a basalt (JMP2), an iron formation (5PLJ009). a volca^suM4PU017) a schist 
(5PLT001) and a granite (5PLJ001) Both weathered (exterior) and unweathered (interior) surface spectra of each sample 
SSrS It Koi ^at physical and chemical weathering operate at low rates in cold, dry teeru The : = 
environment of these regions affects the relative importance of different weathering processes [2). The most co 
alteration phenomona which have been amply documented are the formation of feme iron oxides/hydroxides, as evidenced 
by rust-colored staining on rock surfaces and in soils [3,4.5,6], and the formation of days [3.7.8]. 

Results: Basalt- The weathered and unweathered surface spectra of a basalt arc shown in Figure 1 - The weathered 
surface specmTm shows strong evidence of ferric iron at shorter wavelengths- an « '^orp^ Tre 

broad absorption bands near 0.65 and 0.9„m. The wavelength position of the i inS sX 

consistent with multiple hydrated ferric spedes such as goethite, ferroxyhyte and femhydnte [9,10], The interior surface 
spectrum suggests a greater ferrous iron component and is similar to other unweathered basalt spectra [11.12]. At longer 
w^dSthsK sptSS arc similar. The various absorption bands can be assigned to day lattia-OH absorption bands 
involving both Mg and Al, as expected for a basalt. The presence of these bands in the interior surface spectrum indicates 

that the interior has been slightly metamorphosed. viciHf , w „. hrrjn£ , 

Iron Formation: Sample 5PLJ009 is composed of magneute, chert and sidente. The degree of visible wrather^gis 

small enough so that the underlying banding is still visible. The weathered surface spectrum is 

iron absorption features. The broad band near 1.9ion indicates that weathering has produced multiple 1 hydrated speoes 
and/or proriy -crystallized phases (Figure 2). The interior surface spectrum is dominated by sidente, which exhibits a broad 
absorption band between *1.0 and *l.Vn. However, the expected strong carbonate absorption bands at 2.35 and 2.56nm 

[13] arc absrari^^^ The weathered sur f ace spectrum of a felsic-pyritic metatuff is again dominated by 
at shorter wavelengths (Figure 3). A broad pyrite absorption band expected near lion is not evident The weakness of 
various clay lattice-OH bands suggest that the formation of well-crystalized hydrated species has not been significant. Th 
interior sj^tnimshows almost no evidence Tor ferric iron absorption bands. The broad feature centered near 1.05^ is 

C ^ Ct Tc^° f ThTwMthTr^ a ^d unweathered surface spectra or a cordierite-biotite knotted schist differ “^kedlyfrom 
each other (Figure 4 ). In both cases there is a general lack of diagnostic absorption bands. The general reflectance nse 
towards longer wavelengths and lack of well-defined absorption bands in the extenor spectrum suggests that vrotoemg tas 
generated a number of poorly crystallized phases which absorb incident radiation over a range of energies, probably through 

0,12186 Granite differs from the other samples in its general lack of iron-bearing minerals. This is reflected m the 

lack of ferric iron absorption bands in the spectrum (Figure 5). Both the intenor and extenor surface spectra show th 
cam p types of long wavelength absorption bands, indicative of aluminum. 

Discussion- In spite of the presumed low rate and duration of weathering which has presumably affected these 
samples, exterior (weathered) and interior (unweathered) surfaces of the same sample vary ^^T^S^e^low 
r urall y different (Figures 4 & 5). It is apparent that significant spectral changes can be effected by even the low 
weathering rates prevail in cold desert regions. The formation of ferric iron s^es ^lts tn the app^Mce jf ^ 
absorption edge near 0.53/un, a shoulder or band near 0.64Mm. and an absorption band near O.Vn [9 ^1^. TTiwe fearo^ 
are most apparent in the most iron -rich samples such as basalt, iron formation and volcanoclasuc. The refkcmce spectra 

amply confirm the physico-chemical and spectral importance of ferric iron oxrfe ^^ de J^f ° n J^fiStivelv in 

Aluminum- and magnesium -rich lithologies give nse to Al-OH and Mg -OH absorpuon bands ^respectively in 
weathered surface spectra in cases where the absorption bands are resolvable. For example, the absorption bands present in 
the exterior granite specmim at 2.20, 2.32 and 2.35um can be assigned to Al-OH vibrations. There is no clear whence for 
an Mg-OH absorptionband expected near 2.4*m. This is consistent with the alummum-nch composition of the granite. 

A general lack of well-defined absorption bands attributable to specific hydrated phases indicates that cold desert 
weathering is generally not accompanied by the formation of well -crystallized phases, but that clay formation is nevertheless 
an important process. The broadness of the ferric iron absorption bands also suggests that poorly -crystallized phases are 
important. The formation of poorly -crystallized phases in cold deserts is consistent with interpretations of maruan 
telescopic spectra as indicating semi -amorphous or amorphous phases such as palagomte ■ [14]. , , 

The spectral changes associated with cold desert weathering include variauons in overall spectral slope, the 
appearance/disappearance of absorption bands, shifts in absorption band minima wavelength positions, and changes in band 
shapes and intensities. Dramatic spectral differences between interior and exterior surface spectra are present even when 

weathered surfaces are thin enough so that underlying petrofabrics are still visible. . 

In spite of the small number of samples in the current suite, the data indicate that cold desert weathering o 
carbonates may be sufficient to render them virtually indistinguishable. In addition, the oxidation of ferrous to feme iron 
proceeds regardless of whether the iron is derived from silicates (basalt), oxides (iron formauon), carbonates (iron 
formation), or sulfides (volcanoclastic). It is clear that unweathered, interior rock sample spectra are generally not useful 
for the interpretation of remote sensing data for targets such as Mars which are subject to subaerial weathering. 
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Figure I. Reflectance spectra of interior and exterior surfaces of 
basalt sample JMP2. 



Ftgwe 2. Reflectance spectra of interior and exterior surfaces of 
iron formation sample 5PUO09. ^ 




Flgwe 4. Reflectance spectra of interior and exterior surfaces or 
cordierite-biotite knotted schist sample 5PLT001. 



Figure 5. Reflectance spectra of interior and exterior surfaces of 
granite sample 5PU001. 


Figure 3. Reflectance spectra of interior and exterior surfaces of 
fdsac-pyritic-metatuff (roicanodastic) sample 4PLJ017. 






